The spherical collapse model and cluster formation beyond the A cosmology: 

Indications for a clustered dark energy? 
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We generalize the small scale dynamics of the universe by taking into account models with an 
equation of state which evolves with time, and provide a complete formulation of the cluster viri- 
alization attempting to address the nonlinear regime of structure formation. In the context of the 
current dark energy models, we find that galaxy clusters appear to form at z ~ 1 — 2, in agreement 
with previous studies. Also, we investigate thoroughly the evolution of spherical matter perturba- 
tions, as the latter decouple from the background expansion and start to "turn around" and finally 
collapse. Within this framework, we find that the concentration parameter depends on the choice 
of the considered dark energy (homogeneous or clustered). In particular, if the distribution of the 
dark energy is clustered then we produce more concentrated structures with respect to the homo- 
geneous dark energy. Finally, comparing the predicted concentration parameter with the observed 
concentration parameter, measured for four massive galaxy clusters, we find that the scenario which 
contains a pure homogeneous dark energy is unable to reproduce the data. The situation becomes 
somewhat better in the case of an inhomogeneous (clustered) dark energy. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.+X 



1. INTRODUCTION 

Recent studies in observational cosmology lead to the 
conclusion that the available high quality cosmological 
data (Type la supernovae, CMB, etc.) are well fitted by 
an emerging "standard model". This standard model, 
assuming flatness, is described by the Friedman equation 



a»-l 5 l = 



BttG 



Pm(a) + Px(a)} 



(1.1) 



where a(t) is the scale factor of the universe, p m (a) is the 
density corresponding to the sum of baryonic and cold 
dark matter, with the latter needed to explain cluster- 
ing, and an extra component px{&) with negative pres- 
sure called dark energy needed to explain the observed 
accelerated cosmic expansion (0, 0, H, 0, [E[ and refer- 
ences therein). 

The nature of the dark energy is one of the most 
fundamental and difficult problems in physics and cos- 
mology. Indeed, during the last decade there has been 
a theoretical debate among the cosmologists regarding 
the nature of the exotic "dark energy". Various candi- 
dates have been proposed in the literature, such as a cos- 
mological constant A (vacuum), time- varying A(t) cos- 
mologies, quintessence, k— essence, vector fields, phan- 
toirr tachyons, Chaplygin gas and the list goes on (see 
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and references therein). Within this framework, high en- 
ergy field theories generically indicate that the equation 
of state of such a dark energy is a function of the cosmic 
time. To identify this type of evolution of the equation 
of state, a detailed form of the observed H(a) is required 
which may be obtained by a combination of multiple dark 



energy probes. 

A serious issue here is how (and when) the large scale 
structures and in particular galaxy clusters form. Galaxy 
clusters occupy an eminent position in the structure hier- 
archy, being the most massive virialized systems known 
and therefore they appear to be ideal tools for testing 
theories of structure formation and extracting cosmolog- 
ical information. The cluster distribution basically traces 
scales that have not yet undergone the non-linear phase of 
gravitationally clustering, thus simplifying their connec- 
tions to the initial conditions of cosmic structure forma- 
tion. In the last decade many authors have been involved 
in this kind of studies and have found that the main 
features of the large scale structures (formation epoch, 
shape etc) can potentially affected by the dark energy 
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The aim of this work is along the same lines, attempt- 
ing to investigate the main properties of the non-linear 
spherical model for a large family of dark energy models, 
in which the corresponding equation of state parameter 
varies as a function of the scale factor of the universe, 
w = w(a). The plan of the paper is as follows. In section 
2 we briefly discuss the dark energy issue. In section 3 
we present the various dark energy models and in sec- 
tion 4 we use a joint statistical analysis, in order to place 
constraints on the main cosmological parameters. Sec- 
tion 5 outlines the theoretical analysis of the spherical 
collapse model in which the equation of state parame- 
ter varies with the cosmic time. In section 6 we present 
the corresponding theoretical predictions regarding the 
formation of the galaxy clusters. In section 7, we com- 
pare the predicted concentration parameters with those 
found by four galaxy clusters at relatively large redshifts 
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0.18 < z < 0.45), using the Subaru 8.3 telescope [H|. Fi- 
nally, we draw our conclusions in section 8. Throughout 
the paper we will use H = 70.5km/sec/Mpc. 

2. THE DARK ENERGY EQUATION OF STATE 

In the framework of the general relativity it is well 
known that for homogeneous and isotropic flat cosmolo- 
gies, driven by non-relativistic matter and dark energy 
with equation of state px = w{a)px, the expansion rate 
of the Universe can be written as (see eq ll.ip 



E\a) 



H\a) 



3 d\ny[l+w(y)] 



(2.1) 



Note, that E(a) is the normalized Hubble flow, f2 m is 
the dimcnsionless matter density at the present epoch, 
fix = 1 — O m is the corresponding dark energy density 
parameter and w(a) is the dark energy equation of state. 
Inverting the above equation we simply derive 



w(a) 



-1 



2 Q tflnE 
3 da 



1 - Q m a- 3 E~ 2 (a) 



(2.2) 



The exact nature of the dark energy has yet to be found 
and thus the dark energy equation of state parameter 
includes our ignorance regarding the physical mechanism 
which leads to a late cosmic acceleration. 

On the other hand, it is possible to extent the previous 
methodology in the framework of modified gravity (see 
[13, Ell)- Instead of using the exact Hubble flow through 
a modification of the Friedmann equation we can uti- 
lize a Hubble flow that looks like the nominal one (see 
eq ll.ip . The key point here is to consider that the accel- 
erated expansion of the universe can be attributed to a 
"geometrical" dark energy component. Due to the fact 
that the matter density in the universe (baryonic+dark) 
can not accelerate the cosmic expansion we perform the 
following parametrization [13, [H| : 



E 2 (a) 



H 2 {a) 



a. 



6H Z 



(2.3) 



Obviously, with the aid of the latter approach we include 
any modification to the Friedmann equation of general 
relativity in the last term of eq. (|2.3p . Now from eas. (|2.2[ 
12. 3p we can derive, after some algebra, the "geometrical" 
dark energy equation of state 



w(a) 



1 <M8H 2 
3 dhia 



(2.4) 



From now on, for the modified cosmological models we 
will use the above formulation. 



3. LIKELIHOOD ANALYSIS 

In this work we use various cosmological observations 
in order to constrain the dark energy models, explored 



here (see section 3). First of all, we use the Baryonic 
Acoustic Oscillations (BAOs). BAOs are produced by 
pressure (acoustic) waves in the photon-baryon plasma 
in the early universe, generated by dark matter overden- 
sities. Evidence of this excess was recently found in the 
clustering properties of the luminous SDSS red-galaxies 
[H| and it can provide a "standard ruler" with which we 
can constraint the dark energy models. In particular, we 
use the following estimator: 



A(p) = 



~z 2 s E(a s )]i/3 



da 



a 2 E(a) 



2/3 



(3.1) 



measured from the SDSS data to be A = 0.469 ± 0.017, 
where z s = 0.35 [or a 8 = (1 + z s ) _1 — 0.75] and E(a) = 
H(a)/Ho is the normalized Hubble flow. Therefore, the 
corresponding Xbao function is simply written 



Xbao(p) = 



[A(p) - 0.469] 2 
0.017 2 



(3.2) 



where p is a vector containing the cosmological parame- 
ters that we want to fit. 

On the other hand, a very accurate and deep geomet- 
rical probe of dark energy is the angular scale of the 
sound horizon at the last scattering surface as encoded 
in the location if T of the first peak of the Cosmic Mi- 
crowave Background (CMB) temperature perturbation 
spectrum. This probe is described by the CMB shift pa- 
rameter [45|, 0, H3| which is defined as 




da 



a 2 E{a) 



(3.3) 



The shift parameter measured from the WMAP 5-years 
data Q to be R = 1.71 ± 0.019 at z ls = 1090 [or a ls = 
(l + ^s) -1 ~ 9.17x 10~ 4 ]. In this case, the xlmb function 
is given 



Xcmb(P) 



[i?(p)-1.71] s 
0.019 2 



(3.4) 



Finally, we utilize the Union08 sample of 307 super- 
novae of Kowalski et al. Q. In this case, the XsNia 
function becomes: 



XsNia(p) = 



M th Kp) ~ M° bs (a*) 



(3.5) 



where a,- 



(1 



is the observed scale factor of the 



Universe, Zi is the observed redshift, [i is the distance 
modulus [i = m — M = blogd^ + 25 and dz,(a,p) is the 
luminosity distance 



d L (a,p) = - 

a J a 



V 2 H(y) 



(3.6) 



where c is the speed of light. We can combine the above 
probes by using a joint likelihood analysis: 

Aot(p) = £bAO X ^cmb X £sNIa 
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B. The Braneworld Gravity - BRG model 



xL(p) 



2 2 
XBAO + Xcmb 



XSNIa 



in order to put even further constraints on the parameter 
space used. Note, that we define the likelihood estimator 
1 as: Cj oc exp[~x|/2]. 



4. CONSTRAINTS ON THE FLAT DARK 
ENERGY MODELS 

In this section, we consider a large family of flat dark 
energy models and with the aid of the above observa- 
tional data we attempt to put tight constraints on their 
free parameters. In the following subsections, we briefly 
present these cosmological models which trace differently 
the nature of the dark energy. 



A. Constant equation of state - QP model 

In this case the equation of state is constant (see for 
a review, [15|; hereafter QP-models). Such dark energy 
models do not have much physical motivation. In partic- 
ular, a constant equation of state parameter requires a 
fine tunning of the potential and kinetic energies of the 
real scalar field. Despite the latter problem, these dark 
energy models have been used in the literature due to 
their simplicity. Notice, that dark energy models with 
a canonical kinetic term have — 1 < w < —1/3. On the 
other hand, models of phantom dark energy (w < —1) 
require exotic nature, such as a scalar field with a neg- 
ative kinetic energy. Now using eq. (|2.1[) the normalized 
Hubble parameter becomes 



E 2 {a) = Q r 



(i-a 



-3(l+w) 



(4.1) 



Comparing the QP-models with the observational data 
(we sample Vl m £ [0.1,1] and w £ [—2,-0.4] in steps 
of 0.01) we find that the best fit values are f2 m = 
0.28±0.02andu; = -1.02±0.06 with X tot(^m,w)/dof ~ 
309.2/307 in very good agreement with the 5 years 
WMAP data [4]- Also Davis et al. @ using the 
Essence-SNIa+BAO+CMB and a Bayesian statistics 
found TO = 0.27 ± 0.04, while Kowalski et al. 0] uti- 
lizing the Union08-SNIa+BAO+CMB obtained fl m = 
0.285^q'q2^q oi- Obviously, our results coincide within 
la errors. It is worth noting that the concordance A- 
cosmology can be described by QP models with w strictly 
equal to -1. In this case we find: Q m = 0.28 ± 0.02 with 
xt t (n m )/dof ~ 309.3/308. 



1 Likelihoods are normalized to their maximum values. Note, that 
the step of sampling is 0.01 and the errors of the fitted parameters 
represent 2a uncertainties. Note that the overall number of data 
points used is Ntot = 309 and the degrees of freedom: dof = 
Ntot — "fit, with nflt the number of fitted parameters, which 
vary for the different models. 



In the context of a braneworld cosmology (hereafter 
BRG) the accelerated expansion of the universe can be 
explained by a modification of gravity in which gravity 
itself becomes weak at very large distances (close to the 
Hubble scale) due to the fact that our four dimensional 
brane survives into an extra dimensional manifold (see 
[48| and references therein). The interesting point in this 
scenario is that the corresponding functional form of the 
normalized Hubble flow, eq. (|2.3p . is affected only by one 
free parameter (fl m )- Notice, that the quantity 5H 2 is 
given by 



SH 2 = 2fl h 



(4.2) 



where £l w = (1 — ^m) 2 /4. The geometrical dark energy 
equation of state parameter (see eq |2.4p reduces to 

1 



u>(a) 



1 + £l,„(a) 



(4.3) 



where fl m (a) = £l m a~ J /i? (a). The joint likelihood anal- 
ysis provides a best fit value of Sl m = 0.24 ±0.02, but the 
fit is rather poor Xtoti^m) / 'dof ~ 369/308. 



C. The parametric Dark Energy - CPL model 

In this model we utilize the Chevalier-Polarski-Linder 
( [49|, |50|] , hereafter CPL) parametrization. The dark en- 
ergy equation of state parameter is defined as a first order 
Taylor expansion around the present epoch: 



w(a) — wo + Wi (1 — a) 



(4.4) 



The normalized Hubble parameter is given by (see 
ealO): 



E 2 {a) = n m a-' i + {l-n m )a-^ 1+Wa+Wl) e' Wl(a - 1) (4.5) 

where Wq and W\ are constants. We sample the un- 
known parameters as follows: wq e [—2, —0.4] and w\ £ 
[-2.6,2.6]. We find that for the prior of fi m = 0.28 the 
overall likelihood function peaks at wq = — 1-lln 16 an< ^ 
wi = 0.60^5 54 while the corresponding xfot 
is 307.6/307. 



D. The low Ricci dark energy - LRDE model 

In this modified cosmological model, we use a simple 
parametrization for the Ricci scalar which is based on 
a Taylor expansion around the present time: 1Z = rg + 
n(l — a) [for more details see [43|]]- It is interesting to 
mention that at the early epochs the cosmic evolution 
tends asymptotically to be matter dominated. In this 
framework, we have 



E 2 (a) 



a 



4(r +ri-l) e 4ri(l-a) 



a > a t 
a < a t 



(4.6) 
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where a t — 1 — (1 — 4ro)/4r±. The matter density pa- 
rameter at the present epoch, is directly related with the 
above constants via 



n, 



4rn — 4n 



1 



4n 



4r +4n-l 



(4.7) 



The corresponding equation of state parameter is given 
by 



w(a) 



1-411 r 



1- Q e" Sa^-^idy/y) 



(4.8) 



Notice, that we sample the unknown parameters as fol- 
lows: r G [0.5,1.5] and i\ e [—2.4,-0.1] and here 
are the results: ro = 0.82 ± 0.04 and n = - 
(fim - 0.28) with xto t (r ,ri)/dof ~ 309.8/307 
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E. The high Ricci dark energy - HRDE model 

A different than the previously mentioned geometrical 
method was proposed by Linder [5lj ] . in which the Ricci 
scalar at high redshifts evolves as 



K 



1 — 3wi 



SH 2 



H 2 



(4.9) 



where 5H 2 — E 2 (a) — Q m a> 3 . In this geometrical model 
the normalized Hubble flow becomes: 



E 2 (a) = Q m a- 3 (I + (3a- 3wi ) 



-3«M- lni WMl+/?) 



(4.10) 



where j3 = fi m 1 — 1. Using the same sampling as in the 
QP-models, the joint likelihood peaks at f2 m = 0.28 ± 
0.03 and w x = -1.02 ± 0.1 with Xtoti^m^^/dof ~ 
309.2/307. To this end, the effective equation of state 
parameter is related to E{a) according to eq. 



F. The tension of cosmological magnetic fields - 
TCM model 

Recently, (HJ proposed a novel idea which is based on 
the following consideration (hereafter TCM): if the cos- 
mic magnetic held is generated in sources (such as galaxy 
clusters) whose overall dimensions remain unchanged 
during the expansion of the Universe, the stretching of 
this held by the cosmic expansion generates a tension 
(negative pressure) that could possibly explain a small 
fraction of the dark energy (~ 5 — 10%). In this flat cos- 
mological scenario the normalized Hubble flow becomes: 



E 2 (a) = VL m a 



Qpia 



-3+n 



(4.11) 



where £Ib is the density parameter for the cosmic mag- 
netic fields and = 1 — &m ~~ The equation of 
state parameter which is related to magnetic tension is 
(see ea i2T|) 



w(a) 



30 A + nn B a- 3+n 



fW- 3 +") 



(4.12) 



To this end, we sample Sl^ £ [0,0.3] and n € [0,10] 
and we find that for the prior of fl m = 0.28 the best fit 
values are: £Ib = 0.10 ± 0.10 and n 
X 2 tot (n B: n)/dof ~ 308.9/307. 
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G. The Pseudo-Nambu Goldstone boson - PNGB 
model 

In this cosmological model the dark energy equation of 
state parameter is expressed with the aid the potential 
V(4>) oc [l + cos((f>/F)] ||: 



w(a) 



-1 + (l+w )a l 



(4.13) 



where w G [-2, -0.4] and F S [0, 8]. Obviously, for a < 
1 we get w(a) — > —1. Based on this parametrization 
the normalized Hubble function is given by 



E 2 (a) = n m a~ 3 + (1 - n m ) Px {a) 



(4.14) 



In this context, the corresponding equation of state pa- 
rameter is 



px(a) = exp 



3(1 + w ) 



(1 - a?) 



(4.15) 



Notice, that the likelihood function peaks at wo 
-1.04 ±0.17 and F = 5.9 ± 3.2 with xht(™o, F)/dof 
317/307. 



H. The early dark energy - EDE model 

Another cosmological scenario that we include in our 
paper is the early dark energy model (hereafter EDE). 
The basic assumption here is that at early epochs the 
amount of dark energy is not negligible [55j. In this 
framework, the overall dark energy component is given 
by 



n x ( a ) 



1 - fir, 



n e (i 



-3wn 



n e (i 



-3w \ 



1 ^771 ^771^ 

(4.16) 

where f2 e is the early dark energy density and wo is the 
equation of state parameter at the present epoch. Notice, 
that the EDE model was designed to simultaneously (a) 
mimic the effects of the late dark energy and (b) provide 
a decelerated expansion of the universe at high redshifts. 
The normalized Hubble parameter is written as 



E 2 (a) 



,-3 



l-fijr(a) 



(4.17) 



while using eq. (|2.2p , we can obtain the equation of state 
parameter as a function of the scale factor. Now, from 
the joint likelihood analysis we find that fl e = 0.05 ±0.04 
and w = -1.14+g'^ (for the prior of fl m = 0.28) with 
Xl t (n e ,w )/dof ~ 308.7/307. 
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I. The Variable Chaplygin Gas - VCG model 

Let us consider now a completely different model 
namely the variable Chaplygin gas (hereafter VCG) 
which corresponds to a Born-lnfeld tachyon action [56l . 
[57| . Recently, an interesting family of Chaplygin gas 
models was found to be consistent with the current ob- 
servational data [58j . In the framework of a spatially flat 
geometry, it can be shown that the normalized Hubble 
function takes the following formula: 

E 2 (a) = n b a- 3 + (l-n b )^B s a- 6 + (1 - B s )a~ n (4.18) 

where ~ 0.021/i -2 is the density parameter for 
the baryonic matter [59| and B s e [0.01,0.51] and 
n G [—4,4]. The effective matter density parameter 
is: tt e J f = fi 6 + (1 - Slb)VB~s- We find that the best 
fit parameters are B s = 0.05 ± 0.02 and n = 1.58±o l| 
~ 0.26) with xlt(Bs,n)/dof ~ 314.7/307. 

5. EVOLUTION OF MATTER 
PERTURBATIONS 

In this section we study the spherical collapse model 
by generalizing the basic non-linear equations which gov- 
ern the behavior of the matter perturbations within the 
framework of the previously described dark energy mod- 
els. Also, we compare our predictions with the traditional 
A cosmology. This can help us to understand better the 
theoretical expectations of the current dark energy mod- 
els as well as the variants from the concordance A cos- 
mology. 




FIG. 1: Upper Panel: The deviation (1 - A X /A A )% of the 
collapse factor for various dark energy models with respect 
to the A solution (solid line). Bottom Panel: The density 
contrast at the virialization, A V i r , as a function of redshift. 
The points represent the following cosmological models: (a) 
BRG (open stars), (b) LRDE (open squares), (c) TCM (open 
triangles) , (d) EDE (open circles) and (e) PNGB (solid trian- 
gles). The lines represent: (a) CPL model (long dashed), (b) 
HRDE model (dot line) and VCG model (dashed line). 



A. The Evolution of the linear growth factor 

The evolution equation of the growth factor for models 
where the dark energy fluid has a vanishing anisotropic 
stress and the matter fluid is not coupled to other matter 
species is given by ([42j], [6(3], [6l|, [62j]): 



D 



1 



w(a) 
1 + X(a) 



D 



3 X(a) D 
2 1 + A»a^ 



where 



X(a) 



-3f'w(y)dlny _ ^mg 



1 - a 



(5.1) 



(5.2) 



Note, that the prime denotes derivatives with respect to 
the scale factor. Useful expressions of the grow th factor 
can be found for the ACDM cosmolog y in 16011 . for dark 
energy models with w = const in [63| , [22| , [25| ] , (64| , for 
dark energy models with a time varying equation of state 
in 65] and for the scalar tensor models in [66[. Finally, 
in this work, the growth factor evolution for the current 
cosmological model is derived by solving numerically eq. 
(|5.1|) . Note, that the growth factors are normalized to 
unity at the present time. 



B. The spherical collapse model 

The spherical collapse model, which has a long history 
in cosmology, is a simple but still a fundamental tool for 
understanding how a small spherical patch [with radius 
R(t)] of homogeneous overdensity forms a bound system 
via gravitational instability [67| . From now on, we will 
call at the scale factor of the universe where the over- 
density reaches its maximum expansion (R = 0) and a c 
the scale factor in which the sphere virializes implying 
that a large scale structure has formed, while Rt and R c 
are the corresponding radii of the spherical overdensity. 
Note that in the spherical region, p ms cx R~ 3 is the mat- 
ter density, p m is the background matter density and px s 
denotes the corresponding density of the dark energy. In 
order to address the issue of how the dark energy it- 
self affects the gravitationally bound systems (clusters of 
galaxies), we have to deal with the following situations: 
(i) clustered dark energy, considering that the whole sys- 
tem virializes (both matter and dark energy), (ii) the 
dark energy remains clustered but now only the matter 
virializes and (iii) the dark energy remains homogeneous 
and only the corresponding matter virializes (for more 
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TABLE I: Numerical results. The 1 st column indicates the 
dark energy model used. Between column two and four, we 
present the main properties of the spherical collapse model, 
assuming that galaxy clusters have formed prior to the epoch 
of a c ~ 0.42 (zc ~ 1.4), in which the most distant cluster has 
been found [6g]. Column five, corresponds to the current age 
of the universe, to. Finally, the last two rows correspond to 
the (71, ni) constants which are included in the approximate 
A„, r (a) formula (see eq |5.16[) . 



Model 


Zta 




C 


io/Gyr 


7i 


m 


A 


2.89 


192.0 


5.56 


13.72 








QP 


2.89 


191.3 


5.55 


13.77 


0.384 


2.350 


BRG 


2.98 


247.3 


5.89 


13.40 


0.446 


-0.556 


CPL 


2.90 


204.6 


5.66 


13.62 


0.210 


-1.40 


LRDE 


2.86 


181.2 


5.49 


13.84 


0.152 


-6.848 


HRDE 


2.89 


193.6 


5.63 


13.77 


0.050 


-0.117 


TCM 


2.89 


191.4 


5.56 


13.78 


0.450 


5.387 


PNGB 


2.89 


191.9 


5.57 


13.76 








EDE 


2.88 


200.1 


5.71 


13.71 


0.054 


-1.70 


VCG 


2.85 


185.1 


5.53 


14.01 


0.226 


-1.245 



details see j28|, [23 [3g| and [4(|). Note, that in this 
section we are using the third possibility (for more on 
inhomogeneous dark energy see section 7B). 

In this section, we review only some basic concepts 
of the problem based on the assumption that the dark 
energy component under a scale of galaxy clusters can 
be treated as being homogeneous: pxs(a) — px{&) and 
w s (a) = w(a). In general the evolution of the spherical 
perturbations as the latter decouple from the background 
expansion is given by the Raychaudhuri equation: 



R 
R 



Now, if we perform the transformations 

and y 
then eqs. ()l.ll l5.3|) become: 



a 
a, 



R_ 

Rt 



and 



with 



px.t 

Pm,t 



+ vyl{x,y) 

V 



1 - Q 



rn . / 



m,t 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



where px,t is the dark energy density at the turn around 
epoch and £l m ,t is the corresponding matter density pa- 
rameter. Also, p ms ,t = CPm,t is the matter density for the 
spherical region at the turn around time, while p m ,t de- 
notes the background matter density at the same epoch. 



Note, that in order to obtain the above set of equations 
which govern the global and local dynamics, we have uti- 
lized the following relations: 



Prtis — Prns,t 



Cp 



m,t 



and 



where 



I(x,y)=I{x) = [l + 3w(x)]f(x) 



f( x ) = e 3 ^ dlnu[l+w(u)] 

Finally, fl m (x) is given by 

1 



1 + vx 3 f{x) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



In this context, the time needed for a spherical shell to 
re-collapse is twice the turn-around time, tf ~ 2t t which 
implies 



1 



-da 



1 



-da 



(5.12) 



/q H(a)a J H(a)a 

In the case of the usual A cosmology one can prove that 



sinh 



3/2 



= 2sinh 



3/2 



(5.13) 



where vq = (1 — Q m )/n m . Of course, in order to in- 
clude the dark energy models in our analysis we have 
to solve eq. (|5.12l) numerically. As an example, assuming 
that galaxy clusters have formed prior to the epoch of 
a c ~ 0.42 (z c ~ 1.4), in which the most distant cluster 
has been found [68j, the turn around epoch is not really 
affected by the dark energy component (see Table 1), ic. 
a t ~ 0.26 (or z t ~ 2.8). This is to be expected, due to the 
fact that at large redshifts matter dominates the Hubble 
expansion. It is worth noting that the ratio between the 
scale factors converges to the Einstein de Sitter (fi m = 1) 

value — 2 2 / 3 at high redshifts which implies that 

On the other hand, utilizing both the virial theorem 
and the energy conservation 2 we reach to the following 
condition: 



\r§^(u g + u Xs ) + u G + U Xt 



= [Uc 



U X s} a=a ' 
(5.14) 



2 In fact, a smooth dark energy component violates the energy 
conservation, in the sense that the dark energy freely flows out- 



wards the overdensity I2a, l39ll . Only small scale clustering dark 
energy would satisfy the energy conservation, for instance the 
chameleon models |29l . I30L l3ll . |32| |. Such violation is however 
only important at very late times, when dark energy dominates. 
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where U G = -3GM 2 /5R is the potential energy and 
Uxs — —4:TrGM(l + 3w)px s R 2 /5 is the potential energy 
associated with the dark energy for the spherical over- 
density (see [HI and [IH). Using the above formulation 
we can obtain a cubic equation that relates the ratio be- 
tween the virial R c and the turn-around outer radius Rt 
the collapse factor (A = R c /Rt). Of course in the case 
of w = — 1 the above expressions get the usual form for 
A cosmology ([Hj], [2l[) while for an Einstein-de Sitter 
model (O m = 1) we have A = 1/2. 

Finally solving numerically eq. (|5.14p we calculate the 
collapse factor. In general, we find that the collapse fac- 
tor lies in the range 0.43 < A < 0.50 in agreement with 
previous studies ( [H, [H, 0, [H, [S3, 0). 

Also, in the upper panel of figure 1 we plot the devia- 
tion, (1 — Ax/Aa)%, of the collapse factors Xx(z c ) for 
the current dark energy models with respect to the A 
solution Aa(z c ). It becomes evident that the shape of 
the cosmic structures which produced by the CPL (long 
dashed line), LRDE (open squares), HRDE (dot line), 
TCM (open triangles) and PNGB (solid triangles) mod- 
els is close to that predicted by the usual A cosmology. 
On the other hand, the largest positive deviation of the 
collapse factor occurs for the BRG model (open stars) 
which implies that the latter model produces more bound 
systems with respect to the above dark energy models. 
Therefore, within this cosmological scenario the corre- 
sponding large scale structures should locate in very large 
density environments. The opposite situation holds for 
the VCG (dashed line) and LRDE (open squares) models 
due to their negative deviations. 

In the bottom panel of figure 1 we present the evolution 
of the density contrast at virialization 



pn 



A :i 



(5.15) 



where p ms . c is the matter density in the virialized struc- 
ture while p m , c is the background matter density at the 
same epoch. We verify, that the density contrast de- 
creases with the formation (virialization) redshift z c . Ob- 
viously, the factor A„i r , provided by the spherical col- 
lapse model, plays a key role in this kind of studies. In- 
deed, the differences among the A V i r for the current dark 
energy models enter through the H(z) (see eq. (|5.12p ). 
This feature implies that perhaps the density contrast at 
virialization can be used as a cosmological tool (see sec- 
tion 7). In any case, at very large redshifts the above 
density contrast tends to the Einstein-de Sitter value 
(A V i r ~ 187T 2 ), as it should. In this context, following the 
notations of [27], H^, [7(| , we provide an accurate fitting 
formula to A v i r (within a physical range of cosmological 
parameters) 



where 



187r 2 [l+ 7 0"(a)] 



7 = 0.44-1.31(Ha*)| 71 -1) 



(5.16) 



n = 0.94- 0.21(|w(a*)r 11 - 1) 

a* = 0.4 (or = 1.5) and O(a) = Q^(a)-1. Note, that 
7i and ri\ are constants (see the last rows of Tablel). 

Finally, in Table 1 we list the results of the spheri- 
cal collapse model considering that galaxy clusters have 
formed at z c ~ 1.4, ie., (a) the cosmological models and 
the value of the turn around redshift, (b) the virial den- 
sity Ai,i r , at the collapse time, as well as the density 
excess of the matter density in the spherical overdensity, 
C, at the turn around time, (c) the predicted age of the 
universe and (d) the constants (71, ni) which are included 
in eq. f5TTB|) . 



6. THE FORMATION OF GALAXY CLUSTERS 

In this section we briefly investigate the cluster forma- 
tion processes by generalizing the basic equations which 
govern the behavior of the matter perturbations within 
the framework of the current dark energy models. Also 
we compare our predictions with those found for the tra- 
ditional A cosmology. This can help us to understand 
better the theoretical expectations of the dark energy 
models as well as the variants from the usual A cosmol- 
ogy 
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FIG. 2: The predicted fractional rate of cluster formation 
as a function of redshift for the current cosmological models 
(as = 0.80). 

The concept of estimating the fractional rate of clus- 
ter formation has been brought up by different au- 
thors [7l|, 111]. In particular, these authors introduced 
a methodology which computes the rate at which mass 
joins virialized structures, which grow from small ini- 
tial perturbations in the universe. In particular, the ba- 
sic tool is the Press-Schechter formalism [73j which con- 
siders the fraction of mass in the universe contained in 
gravitationally bound structures (such as galaxy clusters) 
with matter fluctuations greater than a critical value S c , 
[74j . Assuming that the density contrast is normally dis- 
tributed with zero mean and variance a 2 (M, z) we have: 



V(5,z) 



1 



2na(M, z) 



exp 



2cr 2 (M,; 



(6.1) 
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In this kind of studies it is common to parametrize the 
rms mass fluctuation amplitude at 8 ft, Mpc which can 
be expressed as a function of redshift as cr(M, z) = 
ag(z) = D(z)<j%. The current cosmological models are 
normalized by the analysis of the WMAP 5 years data 
as = 0.80 |4j. The integration of eq. ()6.ip provides the 
fraction of the universe, on some specific mass scale, 
that has already collapsed producing cosmic structures 
(galaxy clusters) at redshift z and is given by [72| |: 



P,2, = I 



V(S,z)dS 



P{*) = ~ 2 



1 — erf 



( s c 



\V2a 8 (z) 



(6.3) 



where S c is the linearly extrapolated density threshold 
above which structures collapse, [74|. Notice, that for 
the model of modified gravity (BRG) we use S c ~ 1.47 
(see |zi|), for the EDE model we use 5 C ~ 1.4 (see [76j 
and references therein). To this end, for the rest of the 
dark energy models, due to the fact that w ~ — 1 close to 
the present epoch, we utilize a S c approximation which 
is given by Weinberg & Kamionkowski (see [27J, their 
eq.18). 

Obviously the above generic of form eq. (|6 ,3|) depends 
on the choice of the background cosmology. The next 
step is to normalize the probability to give the number 
of clusters which have already collapsed by the epoch z 
(cumulative distribution) , divided by the number of clus- 
ters which have collapsed at the present epoch [z = 0), 
F(z) — P(z)/P(0). In figure 2 we present in a logarith- 
mic scale the behavior of normalized cluster formation 
rate as a function of redshift for the present dark energy 
models. In general, prior to z ~ 0.2 the cluster formation 
has terminated due to the fact that the matter fluctua- 
tion field, D(z), effectively freezes. For the traditional A 
cosmology we find the known behavior in which galaxy 
clusters appear to be formed at high redshifts z ~ 2 (see 
for example [25[ and references therein) . From figure 2 it 
becomes also clear, that the vast majority of the dark en- 
ergy models seem to have a cluster formation rate which 
is close to that predicted by the usual A cosmology (see 
solid line in figure 2). However, for the BRG and EDE 
cosmological scenarios respectively we find that galaxy 
clusters appear to form earlier (z ~ 3) with respect to 
the CPL, LRDE, HRDE, TCM, PNGB and VCG dark 
energy models. 



COMPARISON WITH OBSERVATIONS 



the profile of cluster dark matter halo which in turn is 
usually modelled by the NFW profile [zl ^ defined as 



Pnfw{R) 



(R/r s )(l + R/r s f 



(7.1) 



where r s and p s are the characteristic radius and density 
of the halo profile. The concentration parameter c V i r is 
an observable characteristic of the halo profile and is de- 
fined as c v i r = r„ir/rs where r V i r is the virial radius of 
the system [zj IZa, UM ■ Dark matter profiles of clusters 
may be determined using lens distortion and magnifica- 
tion [4l], • Such profiles have been recently obtained 
[4l| for four clusters using the wide-field camera Suprime- 
Cam [1U of the Subaru 8.3m telescope. The NFW profile 
fits well the profile of these clusters and the concentra- 
tion parameters have been obtained for all four clusters 

The theoretical prediction of the concentration param- 
eter for each cluster can be made using the results of the 
previous section in the context of the spherical collapse 
model with any dark energy parametrization assuming 
homogeneity of dark energy. In what follows we consider 
the CPL parametrization of equation (|4.4p . In particu- 
lar, we use the approach of Ref. [12] (see also (83|) to 
obtain the predicted value of the concentration parame- 
ter from the density contrast at virialization A„j r . This 
is a simple effective algorithm that approximates well the 
values of the concentration parameter obtained from N- 
body simulations in the context of various cosmological 
dark energy models including the case of the cosmological 
constant. In this approach, the concentration parameter 
c V i r is related to A V i r (z) as 3 



(z c ) n m (« )(i + Zc f 
(z ) n m (z c )(i + z Q ) 3 



(7.2) 



where z c is the formation (virialization) redshift while zq 
is the "observed" redshift (z = zo < z c ). The formation 
redshift z c may be obtained [82j from the relation 



D(z c )a eff (M) 



1 

a, 



(7.3) 



where D(z) is the linear growth factor and a e ff(M) is 
the effective amplitude of the power spectrum on a scale 
M connected with the actual amplitude on a scale M by 



a eff (M) = a(M) 



din a 
rflnM 



(M) 



(7.4) 



C a is a constant with C G ~ 28 [82j . The power spectrum 
amplitude a(M) is well approximated by the fit [H| 



"'-'" = "•1.15 



- 7 i(M)/3 



(7.5) 



A. Homogeneous Dark energy 

A useful observable parameter that is predictable in 
the context of the spherical collapse model is the con- 
centration parameter c V i r . This parameter characterizes 



3 Note that A v i r entering Eq. 117. 2j l is defined with respect to the 
critical density p cr it 18211 , rather than with respect to 

Pm,c as in 

Eq. I|5.15|l . Computing c v i r then requires to multiply A v i r in 
Eqs. 115. 151 1 and (15. 16ft by p m ,o/ Petit- 
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FIG. 3: a: The concentration parameter c V i r (M vir ) for homogeneous dark energy for (wo,Wi) = ( — 1,0) (ACDM) (solid line), 
(wo,Wi) = ( — 1.1,0.6) (best fit of section 4C) and (wo,w\) — (—0.6,0) (dot-dashed line). For comparison we include the 
predicted curve (dashed line) for the values at which the overall likelihood function peaks: (wo,»i) = (—1.1,0.6) (see Sec. [2}. 
Redshift zq = 1.4 has been assumed, b: Same as Fig. 3a with zq = 0. 



where Ms = 6 x 10 14 Qo m /i 1 M® is the mass inside a 
sphere of radius Rg = 8/i _1 Mpc. Also 

7i (M) = (0.3T + 0.2)[2.92 + i log(-^)] (7.6) 

with T = fi 0m /ie- Oi >~ Oi >/° 0m (fi 6 = 0.05 is the baryonic 
density parameter). The linear growth function in equa- 
tion (|7.3p is obtained as (eg [85{ ) 



D(z) 



(7.7) 



where 7 is the effective growth index which depends 
mildly on the dark energy properties in the observation- 
ally allowed parameter range (7 = 0.55 for ACDM). 

Using equations (|7.3|) , (|7.5p , (|7.6[) and the parametriza- 
tion (|5.16p to evaluate A V i r (z) we may use equation (I7.2[) 
to evaluate Ct,i r (M^i r ) for various values of redshift zn and 



compare with the corresponding data of Ref. 41]. The 
results are shown in Fig. 3a for zq = 1.4 and in Fig. 3b 
for zq = 0. In both cases we consider (wq,wi) — (—1,0) 
(ACDM) and {w ,wx) = (-0.6,0) which slightly im- 
proves the fit to the data of Ref. [41[ shown in the same 
figures (values wq < — 1 further decrease the quality of 
fit). We have also tested other cases with 1U1 ^ and we 
have found that the value of c V i r depends (mildly) only 
on the mean value of w(z) at late times when dark energy 
begins to dominate. 

The results shown in Fig. 3 correspond to homoge- 
neous dark energy and confirm the well known puzzle for 
ACDM namely that the predicted concentration param- 
eter in this model in significantly less than the observed 
one [4l[ . The new input of Fig. 3 is that the modification 
of dark energy evolution (within the observationally al- 
lowed range) is unable to resolve this puzzle if homogene- 
ity is assumed for dark energy. As discussed in the next 
subsection however the agreement with observations can 
improve significantly if dark energy is assumed to cluster 
along with the non-linearly clustered dark matter. 



B. Inhomogeneous Dark energy 

Although it has been shown that the dark energy fluid 
does not cluster on scales smaller than lOOMpc |86|, it is 
also interesting to consider the case when the dark energy 
clusters along with the dark matter. Such a behavior is 
expected in models of interacting dark energy with dark 
matter and may also be a result of the gravitational in- 
teraction between nonlinearly clustered dark matter and 
dark energy. In the homogeneous case, the dark en- 
ergy component flows progressively out of the overdensity 
35, 36], and hence energy conservation cannot be applied 
to determine the collapse factor A (along with the virial 
theorem). In order to simplify the inhomogeneous for- 
malism, we consider the extreme case in which the dark 
energy fully clusters along with the dark matter and thus 
we avoid the energy non-conservation problem examined 
in Ref. [35[ . The merit of this assumption is that it al- 
lows an analytical solution to the system of Eqs. (|5.5[) 
and (|5.6p [401 ] due to the fact that the function I(x, y) in 
eq. ([5.6p depends only on y 



I(x,y)=I(y) = [l + 3w(R(y))] 



f(R(y)) 



Therefore, the solution of the system is 



(7. 



Jo 



y 



{ + vyP(y)-({ + P(l)v)y 



1/2 



dy 



[xQ m (x)} 1/2 dx (7.9) 



a t y- 



where P(y) = y 2 / (R(y)) / f \a t ) and R{y) = ( _ 
The matter contrast £ can be obtained by solving nu- 
merically the above equation. 

Although the dark energy clusters, it may not take part 
in the virialization process, in which case the collapse 
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FIG. 4: a: Concentration parameter c vir (M vir ) for the inhomogeneous dark energy when only matter virializes (dashed line) 
and when both matter and quintessence virialize (dot-dashed line) with zo — 0. We include the homogeneous case (solid line) 
for comparison. The three upper curves correspond to (wo,Wi) = (—0.55,0.72), whereas the three bottom ones correspond to 
the values where the overall likelihood function peaks: (100,101) = ( — 1.1,0.60). b: Same as Fig. 4a with zo — 1.4. 



factor is determined by eq. (|5. 14[) . However, if the dark 
energy participates in the virialization, the self-energy of 
the dark energy must also be taken into account. In such 
a case, the equation that determines the collapse factor 
is 135] 



1 dU tot 
2 R -dR +Utot 



= [U tot ] c 



(7.10) 



where 



[/tot, = 



(pm + PXs){<5> m + ^>Xs)dV , 



and (f> m and $x« are the gravitational potentials of the 
matter and dark energy components for a spherical over- 
density. 

After solving for £ and A we use eq. (|5.12[) to compute 
A v i r from eq. (|5.15p at different redshifts, and then c V i r 
from eq. (|7.2p . It was noted in Ref. [36J that the c V i r cor- 
responding to inhomogeneous dark energy models may 
be larger than the corresponding to their homogeneous 
realizations, although such enhancement in c V i r is model 
dependent. Such an increase in c V i r owes primarily to 
the factor A V i r (z c ) / A V i r (zo) in eq. (|7.2p . and it is shown 
in 36] that for certain models this factor may be signif- 
icantly larger than 1. For example, the inhomogeneous 
realization of the model in Ref. [Ill ] results in a c V i r that 
can be larger by a factor of 2. 

In the CPL model we find that, in general, although 
c v i r may be somewhat larger in the inhomogeneous case, 
this improvement is insufficient to account for observa- 
tions if the dark energy component is sub-dominant at 
the time of turn-around. The situation improves when 
we allow the dark energy to give a non-negligible contri- 
bution at turn-around, which can be arranged by taking 
w i w o (see FigHJ). As mentioned before, in such a 
case energy conservation cannot be used to determine 
the collapse factor. However, it is interesting to examine 



the predictions for c vir obtained by pretending that en- 
ergy is conserved, hence using eqs. (|5.14[) and (|7.10[) to 
determine A. Of course, in this case one must understand 
c v i r as an estimate to its actual value. 

Assuming then energy conservation to estimate A, we 
begin by the case when only the matter component viri- 
alizes. In that case A V i r (z c ) / A„ r (zo) in eq. (|7.2p is of 
order 1, and so the enhancement factor in c V i r is lost. 
We find that in the range of M V i r of interest to observa- 
tions: 1O 15 M < M vir < 10 16 M Q [4l|, the concentration 
parameter behaves roughly as c V i r ~ T+f £ ' The agree- 
ment with observations then depends almost exclusively 
on z and z c , requiring an early rapid collapse of the halo 
subunits: z c ~ 7 for zq ~ 0. The agreement with observa- 
tions ameliorates if the dark energy component also viri- 
alizes. The reason is that in such a case the correspond- 
ing collapse factor is smaller than if only matter virial- 
izes [13]. Consequently, the factor A vir (z c )/ A vir (z ) in- 
creases, becoming significantly larger than 1 in the obser- 
vational range of M V i r provided Wi is properly tuned (see 
Fig. H} . Similar results to the ones presented in Fig. 0] 
may be obtained for a wide range of values of wo (as- 
suming of course w± is conveniently tuned) , in particular 
for Wq ~ — 1, which is in better agreement with current 
observations. In FigUJb we present the results for c V i r 
for zq — 1.4. In this case, the factor A V i r (z c ) / A V i r (zo) 
is counteracted by the term \^_ z z a , and hence agreement 
with observations cannot be attained in this case. 

To conclude we would like to emphasize that the re- 
sults presented in Fig. 31 owing to the aforementioned 
energy non-conservation problem, must be understood 
as an estimate to the actual value of c V i r - A more re- 
liable estimate for c V i r must obviously account for en- 
ergy non-conservation in the equation determining the 
collapse factor A. As pointed out in [3^], taking this 
effect into account results in a further reduction of the 
collapse factor A. Therefore, one can expect that a more 
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detailed study will place c V i r in better agreement with 
observations. 



8. CONCLUSIONS 

In this paper we study analytically and numerically the 
spherical collapse model beyond the usual A cosmology, 
by using several parameterizations for the dark energy 
equation of state. In this framework, we first perform a 
joint likelihood analysis in order to put tight constraints 
on the main cosmological parameters by using the cur- 
rent observational data (SNIa, CMB shift parameter and 
BAOs). For the vast majority of the dark energy mod- 
els, we find that the large scale structures (such as galaxy 
clusters) start to form prior to z ~ 1 — 2 [13, [zH ■ The 
amplitude and the shape of the concentration parameter 
(c vir ) is affected by the considered status of the dark en- 
ergy model (homogeneous or clustered). We verify that 
in the case where the distribution of the dark energy 
is clustered we are producing more concentrated struc- 



tures with respect to the homogeneous dark energy. Fi- 
nally, we go a step further by comparing the predicted 
concentration parameters with those observed for four 
massive galaxy clusters (0.18 < z < 0.45) and we find 
that the homogeneous dark energy pattern is unable to 
reproduce the data. The situation becomes somewhat 
better in the case of an inhomogeneous (clustered) dark 
energy. In order to confirm such a possibility, we need to 
create robust cluster surveys at relatively high redshifts 
(1.5 > z > 0.2). The latter result points to the direc- 
tion that perhaps the c V i r parameter can be used in the 
future in order to understand the physical properties of 
the dark energy. 
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